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quantitative  Schlieren  system  used  to  measure  the  radial  density  profiles  of  reduced  density 
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quire  a  coherent  light  source.  The  cylindrical  symmetry  of  the  channels  permits  an  Abd-tikc  Inver¬ 
sion  of  the  path  integrated  bending  of  the  tight  to  yield  the  radial  density  profile.  The  inversion 
ilgorithm  and  a  program  using  It  are  given  and  some  experimental  results  are  shown, 
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A  QUANTITATIVE  SCHLIEREN  SYSTEM  FOR  MEASURING  RADIAL 
DENSITY  PROFILES  OF  REDUCED  DENSITY  CHANNELS  IN  GASES 

I.  INTRODUCTION 

Hot,  reduced-density  channels  through  gaseous  atmospheres  are  of  interest  for  a  variety  of  rea¬ 
sons.  In  future  inertial  confinement  fusion  reactors  such  channels  will  probably  be  used  to  transport 
the  charged  particle  beams  through  the  reactor  atmosphere  to  the  pellet.1  2  While  in  the  atmosphere 
every  flash  of  lightning  creates  a  hot,  reduced-density  channel  and  it  is  in  the  details  of  the  process  of 
channel  cooling  that  atmospheric  nitrogen  is  fixed.3  In  experiments  at  NRL4,5  pertinent  to  both  these 
areas,  the  channels  produced  in  the  atmosphere  by  laset -guided  and  unguided  electric  discharges  have 
been  studied  using  a  quantitative  Schlieren  technique.6 

This  technique  detects  the  total  angular  deflection  of  light  passing  through  the  channel  as  a  func¬ 
tion  of  the  lateral  distance  from  the  channel  axis.  An  Abel-like  inversion  is  applied  to  find  the  radial 
density  gradient  as  a  function  of  radius  within  the  channel.  The  density  gradient  is  then  integrated  to 
yield  the  density  profile.  The  ability  of  this  system  to  detect  angular  deflections  is  dependant  on  a 
geometric  ray  approximation  being  valid  throughout  the  optical  train.  This  implies  a  limit  on  the  nar¬ 
rowness  of  the  slit  used  in  the  system  which  in  turn  implies  a  minimum  readable  angular  deflection  of 
-  10"4  radians.  This  sensitivity  is  adequate  to  measure  the  strong  density  gradients  encountered  in  the 
channels.  On  the  other  hand,  the  total  deflections  encountered  are  sufficiently  small  that  the  light  rays 
may  be  approximated  as  passing  straight  across  the  channel  when  the  incremental  deflections  along  the 
path  are  summed. 
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II.  A  DESCRIPTION  OF  THE  SYSTEM 

The  quantitative  Schlieren  system  is  shown  in  Fig.  I.  The  ruby  laser  produces  a  parallel  beam  of 
light  ~  1  cm  in  diameter.  The  laser  energy  is  ~  2J  and  the  pulse  width  is  ~  25  ns  FWHM.  The  first 
cylindrical  lens  (a)  brings  the  beam  to  a  horizontal  line  focus  (b).  The  light  then  re-expands  vertically 
until  it  impinges  on  the  first  large  lens  (c)  as  a  slender  vertical  elipse.  This  lens  is  positioned  at  its  focal 
distance  from  the  line  focus  (b)  and  stops  the  vertical  expansion  of  the  light.  At  the  same  time  it 
causes  the  light  to  focus  sideways  to  produce  a  vertical  line  focus  (d).  The  light  re-expands  horizontally 
beyond  the  vertical  line  focus  until  falling  as  another  slender  vertical  elipse  on  the  symmetrically  placed 
second  large  lens  (e).  This  lens  restores  the  horizontal  parallelism  of  the  light  and  focuses  it  vertically 
to  produce  a  second  horizontal  line  focus  at  the  inclined  slit  (j).  The  slit  has  a  width  of  —.5  mm  and 
light  rays  preserve  their  identity  in  transiting  the  slit.  The  central  portion  of  the  horizontal  line  focus 
emerges  from  the  slit,  expanding  vertically,  and  passes  through  the  center  of  the  second  cylindrical  lens 
(g).  This  lens  does  not  affect  the  vertical  expansion  but  projects  an  enlarged  version  of  the  horizontal 
extent  of  the  emergent  light  on  to  the  film  (h).  This  results  in  a  vertical  band  of  light  on  the  film  (Fig. 
2). 


Wnen  a  disturbance  is  present  (/),  the  light  at  various  offsets  may  be  bent  vertically:  such  rays  are 
indicated  by  dotted  lines  in  Fig.  1.  These  rays,  because  of  their  angular  deviation,  reach  a  horizontal 
line  focus  at  the  inclined  slit  (J)  which  is  above  that  formed  by  the  undeviated  rays.  Therefore,  a  non¬ 
central  portion  of  this  light  emerges  from  the  slit  and  is  projected  sideways  by  the  second  cylindrical 
lens.  The  vertical  path  of  these  rays  is  nearly  the  same  beyond  the  slit  as  if  they  were  not  deflected  so 
they  fall  in  their  appropriate  position  vertically  on  the  film.  The  image  on  the  film  becomes  therefore  a 
plot  of  angular  deflection  versus  offset.  For  a  cy lindrically  symmetric  channel  there  will  result  a  record 
which  is  antisymmetric  about  the  channel  center  (Fig.  2). 

The  magnification  of  the  offset  distance  depends  on  the  ratio  of  the  focal  length  (/>)  of  the 
second  large  lens  <e0  and  the  distance  (d i)  from  the  inclined  slit  (/)  to  the  film  (//) . 


2 


v  jAa-A-  * 


NRL  MEMORANDUM  REPORT  439C 
■Pfilm  <h 
-Pchannsl  fl 


(1) 


The  angular  magnification  is  determined  by  the  focal  length  (/>)  of  the  second  large  lens  (e),  the 
angle  of  inclination  (<f>)  of  the  slit  (J),  the  focal  length  (/3)  of  the  second  cylindrical  lens  (g),  and  the 
distance  ( d-> )  from  the  second  cylindrical  lens  (g)  to  the  film  (h). 


xfilm 
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cm 
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The  offset  magnification  may  be  experimentally  calibrated  by  inserting  an  opaque  object  of  known 
size.  The  angular  magnification  may  be  experimentally  calibrated  by  inserting  either  a  glass  wedge  of 
known  deviation  or  an  extremely  long  focal  length  lens. 


III.  DATA  ANALYSIS 


The  rate  of  deflection  of  a  light  ray  traveling  in  the  z  direction  due  to  a  transverse  gradient  in  the 
refractive  index  idy/By)  is  given  by7 


(3) 


<10  ^  J_  8t i 
dz  t)  By' 

In  the  channels  being  studied  at  NRL  the  refractive  index  is  determined  by  the  density  cf  polarizable 


neutral  molecules  and  atoms  and  is  near  unity.  In  that  case  it  may  be  shown  that  the  Claussius- 
Mossotti  relationship  between  the  refractive  index  (77)  and  the  species  polarizability  («)  reduces  to8 


7} —  l  ~  knia  (4) 

where  n  is  the  species  number  density  -  We  may  therefore  relate  the  rate  of  deflection  to  the  transverse 


gradient  in  number  density. 
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Valuer  of  2:ra  are  given  in  Table  l  for  several  gases.9 10  n  u  The  polarizabilities  of  a  molecule  and  an 


atom  for  normally  diatomic  gases  obey 
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For  the  situations  considered  here,  the  fractional  dissociation  (8)  is  <  0.2  and  for  simplicity  we  have 
taken 

a{M)  —  a  (A/2). 

which  results  in  an  underestimate  of  the  density  on  the  channel  axis  by  ~  10%. 

We  will  assume  the  reduced  density  channe.s  are  cylindrically  symmetric  and  that  the  light 
traverses  them  along  straight  chords.  We  will  also  assume  the  channel  consists  of  concentric  shells, 
each  characterized  by  a  constant  radial  density  gradient  (Fig.  3).  In  this  diagram  each  shell  is  identified 
by  its  inner  radius  so  that  (bn/br),  represents  the  constant  radial  density  gradient  characterizing  the 
shell  with  inner  radius  rt.  The  outside  of  the  channel  has  radius  rv  and  the  innermost  cylinder  is 
assumed  to  be  uniform  giving  N-  1  nontrivial  shells.  A  ray  is  identified  by  the  radius  of  the  interior 
shell  to  which  it  is  tangent  (rj).  The  radius  is  equal  to  the  offset  of  the  ray. 


The  incremental  deflection  along  a  chord  within  one  shell  may  be  evaluated  analytically  (Fig.  4). 
Because  the  channel  is  cylindrically  symmetric,  the  total  benorng  of  the  light  ray  when  passing  both  into 
and  out  of  a  shell  may  be  found  by  integrating  bO/b:  along  the  inbound  chord  segment  and  doubling 
the  result,  It  is  convenient  to  use  r  as  the  variable  of  integration.  The  variables  /  and  r  are  related  via 

/  =  (r’_  ,/)*  -  -  /•/)*  (6) 

and 

dl  =  (r2  -  r  dr.  (7) 

The  transverse  density  gradient  is  related  to  the  radial  gradient  via 


bn 
0  y 


sin  >b 


bn 
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The  path  integral  for  the  total  bending  in  one  shell  is  therefore  given  by 


10  =  2  J*  I?-  (//  =  2  X 


2 rra  r,  ft,, 

r  dr 

r  br 

ylr-r; 

(8) 


NRL  MEMORANDUM  REPORT  4390 


=  4t ra  r-, 


Bn  Cri- 1  dr 


,  Bn  ,  0+i  +  Vfr+i  “  0 
477a  r,  —  In  - 7=== — 

3rJ/  o  +  V  r,2-/-/ 


We  may  evaluate  the  total  deflection  of  the  ray  tangent  to  r}  (see  Fig.  3)  by  summing  the  bending 
in  each  shell  traversed  by  the  ray 

,>  .<  Bn  i  0+ 1 4*  "v/o+i  —  0”  / 1 

=  477a  rj  £  —  In - j= - —  .  (10) 

i-j  r,  +  VO  “  0 


l |  n  +  y/rf-rf  j 

We  may  find  an  inversio'  of  this  equation  by  dividing  by  r}  and  separating  the  term  representing  the 
innermost  contribution  to  the  bending 


n  ,  f  0+i  +  V 0*i  “02 1  .  ,  ,  O+i  +  V o+i  *• 


0+i  + a/ O+i  -0 


o+Vo2“0" 


1 

iU/v  r  In 

o+i+Vo+i  ~  ri 

O+i+Vo+i  ~r} 

1 

} 

wa  fj  in 

0 

in 

0 

This  equation  allows  us  to  find  ( Bn/Br)j  given  0,  provided  we  already  know  ( Bn/Br ),  for  all 
t  >  j.  We  may  therefore  find  all  values  of  (Bn/Br)  by  successive  solutions  of  this  equation  starting 
with  j  =  .V- 1  and  working  inwards.  This  inversion  forms  the  basis  of  the  computer  program  used  to 
analyze  the  data.  Having  found  Bn/Br  as  a  function  of  radius,  integration  yields  the  density  as  a  func¬ 
tion  of  radius. 

n(r)  =  «(«)  -  f  ~  dr.  (13) 

J  r  Br 

The  computer  program  carrier  this  integration  inwards  one  shell  at  a  time  using  the  assumption  that 


(Bn/Br)  is  constant  within  each  shell. 
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A  listing  of  the  computer  program  is  given  in  ihe  Appendix. 

By  applying  equation  (9)  to  the  piecewise  analytic  density  profile  shown  in  Fig.  5,  a  set  of  test  data 
was  generated  to  check  the  inversion  program.  The  fit  achieved  by  the  program  when  run  on  the  test 
data  is  indicated  on  the  figure. 

A  refinement  which  has  been  added  to  the  inversion  program  is  a  correction  for  the  small  lateral 
position  error  that  results  when  the  light  is  deflected.  This  aberration  is  seen  in  Fig.  2  where  the 
deflected  line  overhangs  the  undisturbed  portions.  The  size  of  this  aberration  is  proportional  to  the 
angular  deflection  and  the  raw  data  is  thus  readily  corrected  prior  to  inversion.  The  correction  to  the 
lateral  position  is  found  as  follows 

■Vreal  *="  ^measured  "b  6  (y)  J  2  (15) 

where  bending  away  from  the  channel  represents  a  positive  angle  and  is  the  type  of  deflection  given  by 

a  reduced  density  cylinder. 


A  test  which  has  been  made  on  the  density  profiles  from  the  ohmicly  heated  channels  is  to  count 
the  total  number  of  equivalent  molecules  in  a  cylindrical  volume  which  contains  the  disturbed  region 
and  to  compare  this  to  the  number  of  molecules  tn  an  equal  size  cylinder  of  ambient  air.  Account  must 
be  taken  of  the  dissociation  in  the  hot  inner  parts  of  the  channel.  The  fractional  dissociation  (S)  is 
defined  by  the  relation 

{!+£)-=  —  (16) 
m 

where  n  is  the  particle  density  and  m  is  the  molecular  density.  A  plot  of  (1  +8)  versus  T  for  air  in 
equilibrium  at  elevated  temperatures  is  shown  n  Fig.  6.  An  approximation  to  this  data  is  shown  by  the 
dotted  line  and  is  given  by 


*9.362  x  10"5  T 


(l  ri)  « 


1 


.751  T  >  2650°* 
or 

T  ^  2650  °A 


(17) 


We  are  interested  in  density  profiles  ut  times  when  pressure  equilibrium  is  established  across  the  chan¬ 
nel.  Therefore  within  the  channel  we  may  write 


nkT  =  1  atm. 


(18; 
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Equations  16,  17,  and  18  represent  three  equations  for  the  unknowns  m,  (1  +  8),  and  T.  We  may 
eliminate  (1+8)  and  T  resulting  in  the  following  solution  for  m. 


m 


the  lesser  of 
n2 


6.875  x  1017  +  .751  n 
or 
n 


(19) 


Using  this  equation  we  convert  the  particle  density  («)  .r;  a  function  of  radius,  into  an  equivalent 
molecule  density  On).  In  the  over  pressured  region  behind  the  shock  wave  application  of  Equation  19 
assigns  zero  dissociat  on.  We  then  count  up  the  equivalent  molecules  per  unit  length 

M  =  JQ '  m  (r)  hr  r  dr.  (20) 

The  program  carries  this  integiation  inwards  one  shell  at  a  time  assuming  m{r)  varies  linearly  within 


each  shell. 


Mj  »  M  *i  +  -y-£(/)+|  -  r?)  +  rrh  (r£.|  +  rf) 

where  g  is  the  slope  of  the  line  *-  variation  in  m(r)  and  li  is  the  intercept. 

„  'Vi  ~  m> 
o+i-o 

li  “  n\j  -  gr,. 

The  program  also  computes  the  comparison  value  of  rrr$m0  where  m0  is  the  ambient  value. 


(21) 

(22) 

(23) 


IV.  SOME  EXPERIMENTAL  RESULTS 


Figure  7  shows  the  density  profile  that  results  when  the  data  given  in  Fig.  2  is  inverted.  The 
figure  also  shows  the  result  of  a  hydrocode  simulation11  for  two  possible  initial  conditions.  The  actual 
>.hun.nel  is  not  perfectly  symmetric  as  ma,  be  seen  from  the  irregular  outline  in  the  Schlieren  photo¬ 
graph  u  Fig.  2.  The  channel  also  becomes  turbulent  with  time  which  indicates  the  presence  of  vorticity 
within  .he  channel.  This  lack  of  symmetry  results  in  errors  which  accumulate  as  the  inversion  proceeds 
inward,  f  le  system  cannot  therefore  resolve  the  low  M/20  atmospheric)  interior  density  predicted  by 
theory.  However  shape  is  predicted  correctly  down  to  densities  of  —  1/10  atmosphere  and  the  overall 
agreenu-.i.  oetween  calculated  and  measured  aensity  profiles  is  well  within  the  uncertainly  in  the  input 
data  (energy  deposited,  deposition  profile,  and  time)  fed  to  the  computer  simulation. 
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V.  CONCLUSIONS 

The  quantitative  Schlieren  technique  described  in  this  paper  is  best  adapted  to  measuring  the  den¬ 
sity  profiles  of  highly  symmetric  (cylindrical  or  spherical)  disturbances  in  gases.  It  is  a  low  sensitivity 
system  useful  for  viewing  disturbances  with  strong  density  gradients.  The  ability  to  resolve  low  abso¬ 
lute  interior  densities  is  dependent  upon  the  degree  cf  symmetry  and  lack  of  interior  structure  (e.g., 
shock  waves  or  turbulence)  of  the  disturbance.  Disturbances  likely  to  meet  these  requirements  include 
electric  sparks,  laser  breakdowns,  and  channels  formed  by  intense,  relativistic  particle  beams,  particu¬ 
larly  those  possessing  a  smooth  (i.e.,  Bennett)  profile. 

The  technique  directly  compliments  normal  Schlieren  pnotography.  It  identifies  the  sign  and  mag¬ 
nitude  of  the  angular  light  deflections  leading  to  darkened  regions  in  a  normal  Schlieren  photograph. 
This  may  be  used  to  identify  areas  of  compression  and  rarefaction  without  actually  inverting  the  data. 
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VII.  APPENDIX:  A  Basic  Computer  Program  for  Inverting  Schlieren  Data 
in  this  program: 

.4  S  is  simply  a  name  for  data  identification 

\1\  is  the  offset  magnification  as  defined  by  Eq.  (1). 

Ml  is  the  angular  magnification  as  defined  in  Eq.  (2). 

VI  is  the  number  of,  pairs  of  data  points,  R(J)  and  0(J ),  where  R(J),  is  'he  offset  on  the 
record,  yfi|m,  and  OU)  is  the  deflection.  Xfiim-  These  are  shown  in  figure  5b. 
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Test  data  is  included  with  the  program  for  the  analytic  model  in  figure  5,  we  note  that  for  this  modtl 
the  value  of  tt r$  m0  is  meaningless  and  should  not  be  compared  with  the  numbers  of  equivalent 
molecules  within  the  raduis  rv. 


80/10/01.  11.44,52. 

PROGRAM  JRGMR 

00100  REM  CYLINDRICAL  LENS  SCHLIEREN  INVERSION 
00110  REM  INPUT  A$,M1,M2 

00120  REM  A$=RECORD  NAME,M1=RADIUS  MAG,M2=ANGLE  MAG  (CM/RADIAN) 
00130  REM  DATA  INPUT  N1,R(J),0(J) 

00140  REM  N1=NO  PTS,R(J)=RAD  (CM) ,0(J)=ANGLE  (A; 

^0150  REM  DATA  INPUT  WITH  INCREASING  RADIUS 
'0160  DIM  R (100) ,0(100) ,X( 100) ,D( 100) ,N( 100) ,P( 100) 

00170  READ  A$, Ml, M2, Ml 
00180  FRINT  "FOR  RECORD”  A$ 

00190  FOR  1=1  TO  N1 
0020?  READ  R(I) ,  0(1) 

00210  0(I)=0(I)/M2 

00220  R(I)=(R(I)/M1)+76. 2*0(1) 

00230  REM  AB3ERATI0N  CORRECTION  INCLUDED 

00240  NEXT  I 

00250  R(0)=0 

00260  R(N1+1)=R(M1)*2 

00270  N(N1)=2.55E19 

00280  P(N1 )=0 

00290  M=0 

00300  R1=R(N1 ) 

00310  R2=R(N1+1 ) 

00320  X(M1  )=0(N1  )/(2*R1*L0G(  (R2+SQR(R2''2-Rr2)  )/R1 ) ) 

00330  FOR  JsM1-1  TO  2  STEP  -1 
00340  R1=R(J) 

00350  R4=R(J+1 ) 

00360  Y=0 

00370  C=LCG((R4+SQR(R4,'2-Rr2))/R1) 

00380  FOR  IsMI  TO  J+1  STEP  -1 
003S0  R3=R(I+1) 

00400  R2=R(I) 

00410  Y=Y+X(I)*LCG((R3+SQR(R3''2-Rr2))/(R2+SQR(R2"2-Rr2))) 
00420  NEXT  I 

00430  X(J)=0(J)/(2#R1*C)-Y/C 

00440  D ( J ) =X( J ) /1 0 . 84E-24 

00450  REM  INTEGRATE  DENSITY  DERIVITIVE 

00460  N(J)=N(J+1)-D(J)*(R2-R1) 

00470  REM  EQU.  MOL.  COUNT 

00480  P(J)=N(jr2/(.751*N(J)+6.875E17) 

00490  IF  P(J)<N(J)  THEN  00510 
00500  P(J)=N(J) 

00510  G=(P(J+1)-P(J))/(R2-R1) 
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00520  H=P(J )~G*R1 

00530  M=M+2 . 094* (R2*3-R 1 "3 ) *G+3 . 1 41 6* (R2~2-R r 2 ) *H 
00540  NEXT  J 
00550  N(1)=N(2) 

00560  N(1 )=N(2) 

00570  M=K+3 . 1 4 1 6*P  ( 2 ) *R  1  "'2 
00580  REM  COMPARISON  VALUE 
00590  0=2.55E19*3.1416*R(Nir2 
C 0600  PRINT  USING  00610 
00610  : RADIUS  DENSITY 

00620  :#.## .  1LM . 

00630  :CM  CM-3 

00640  FOR  J=1  TO  Ml 

00650  PRINT  USING  00620  ,R(J),N(J) 

00660  NEXT  J 
0067 C  PRINT  USING  00630 
00680  PRINT  USING  00700 
00690  PRINT  USING  00710  ,M,Q 
00700  :EQU.  MOLS.  COMPARE 

00710  j//,////—  //.//// . 

00720  REM  DATA  IMPUT 

00730  DATA  "  ANALYTIC  ",2.5,2800,21 

00740  DATA  0, 0,. 243,. 103,. 486,. 208,. 728,. 317 

C0750  DATA  .971,. 432, 1.212,. 555, 1.453,. 692, 1.692,. 849 

00760  DATA  1.929,1.040,2.162,1.298,2.367,1.959,2.616,1.971 

00770  DATA  2.867,1.961,3.119,1.926,3.373,1.865,3.629,1.774 

°0780  DATA  3. 888, 1.649, 4. 149, 1.480, 4. 415, 1.250, 4. 688,. 912 

00790  DATA  5.000,0 

00800  END 

READY. 

RUN 


80/10/01. 

11.46.22. 

PROGRAM 

JRGMR 

FOR  RECORD 

ANALYTIC 

RADIUS 

DENSITY 

.00E+000 

1.01E+018 

1.00E-001 

1.01E+018 

2.00E-001 

1 .00E+018 

3.00E-001 

9.S9E+017 

4.00E-001 

9.99E+017 

5.00E-001 

9.99E+017 

6.00E-001 

9.99E+017 

7.00E-001 

9.99E+017 

3.00E-001 

9.99E+017 

9.0CE-001 

1 .00E+018 

1 . 00E+000 

1.00E+018 

1.10E+000 

3.45E+018 

1 . 20E+000 

5.90E+018 

1.30E+00C 

8.35E+018 

1.40E+000 

1.08E+019 
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1 . 50E+000 
1.60E+000 
1.70E+000 
1.80E+000 
1.90E+000 
2.00E+000 
CM 

EQU.  MOLS. 
1 . 24E+020 


1.33E+019 

1.57E+019 

1.82E+019 

2.06E+019 

2.31E+019 

2.55E+019 

CM-3 

COMPARE 

3.20E+020 


SBU  0.732  'JNTS. 
RUN  COMPLETE. 
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Table  1  —  Polarizabilities  of 
Various  Species 


Species 

2tt  a  (cm2) 

Ref. 

h2 

5.14  x  10~ -4 

9 

H 

6.62  x  10"24* 

10,11 

NS 

1.104  x  10'22 

9 

N 

7.1  x  10" 24 

12 

O, 

1.06  x  10- 22 

9 

0 

4.8  x  10~24 

12 

Air 

10.84  x  10~24 

9 

'estimated  based  on  oscillator  frequencies 
&  strengths 
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Tig  2  -  Data  for  a  laser  initiated  ohmicly  heated  reduced  density  channel,  a)  A  Schliercn  photograph  taken  with  a  pinhole  in 
place  of  the  usual  knife  edge  b)  The  quantitative  Schlieren  photograph  showing  deflections  due  to  both  the  reduced  density 
channel  and  the  outward  going  shock  wave  Both  photographs  were  taken  —  85ycs  after  the  high  voltage  discharge  which  causes 
the  ohmic  heating.  Exposure  time  for  both  photographs  was  —25  ns. 
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